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ABSTRACT 

Ongoing studies with XMM-Newton have shown that powerful accretion disc winds, as re¬ 
vealed through highly-ionised FeK-shell absorption at £ ^ 6.7 keV, are present in a signih- 
cant fraction of Active Galactic Nuclei (AGN) in the local Universe ( [Tombesi et al.|2010a| l. In 
Gofford et al. ( 2013| l we analysed a sample of 51 Suzaku-ohserved AGN and independently 
detected Fe K absorption in ~ 40% of the sample, and we measured the properties of the ab¬ 
sorbing gas. In this work we build upon these results to consider the properties of the associ¬ 
ated wind. On average, the fast winds (vw > 0.01 c) are located (r) ^ 10*^^'^ cm (typically ^ 
10^^^ Ts) from their black hole, their mass outflow rates are of the order (M, 
yr^* or ^ (0.01 — l)MEdd and kinetic power is constrained to {L^. 
lent to ^ (0.1 — 10%)LEdd- We find a fundamental correlation between the source bolometric 
luminosity and the wind velocity, with Vw and a — 0.4+q 2 (90% conhdence), which 

indicates that more luminous AGN tend to harbour faster Fe K winds. The mass outflow rate 
Mw, kinetic power Lw and momentum flux />w of the winds are also consequently correlated 
with Lboh such that more massive and more energetic winds are present in more luminous 
AGN. We investigate these properties in the framework of a continuum-driven wind, show¬ 
ing that the observed relationships are broadly consistent with a wind being accelerated by 
continuum-scattering. We And that, globally, a significant fraction (~ 85%) of the sample can 
plausibly exceed the Lw/Tfioi ~ 0.5% threshold thought necessary for feedback, while 45% 
may also exceed the less conservative ^ 5% of Lboflhreshold as well. This suggests that the 
winds may be energetically significant for AGN-host-galaxy feedback processes. 
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1 INTRODUCTION 

Outflows of photo-ionised gas are ubiquitous in Active Galactic 
Nuclei (AGN). Recent surveys with XMM-Newton and Suzaku by 
ITombesi et al.| ( |2010bt and [Gofford et al.| ( [2013^ (hereafter ‘Pa¬ 

per F) have shown that FeK-shell outflows are observed in a 
significant fraction (~ 40 — 50%) of active galaxies. These po¬ 
tentially massive outflows (Ah ~ 10^^ cm^^), which are identi¬ 
fied through blue-shifted resonant absorption from FeXXV Hea 
(£rest = 6.7 keV) and Fe XXVI Lya (£re.st = 6.97 keV), can posses 
outflow velocities (Vw) ranging from a few thousand kms^* up 
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to mildly relativistic values (Vw ~ 0.1 — 0.3 c; [P ounds et al.|2003[ 
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In Paper I we performed a systematic search for Fe K absorp¬ 
tion in 51 ^Mzafen-observed AGN which were heterogeneously se¬ 
lected from the HEASARC data archiv^(see Paper I for selection 
criteria). Our main results from that paper were: (i) blue-shifted 
Fe K-shell absorption is present at the ^ 2 ct level in 20/51 (~ 40%) 
of the selected AGN (in 28/73 individual observations); (ii) the ab¬ 
sorbing gas is typified by mean values of log(AH/cm^^) ~ 23 and 
log(i^/ergcms^^) ~ 4.5, while the outflow velocity spans a con¬ 
tinuous range between v* < 1,500 km s^ * and < 100,000 km s^ *, 
respectively, with a mean and median velocity of 0.1 c and 0.056 c. 
Ultimately, these results are all consistent with those found by 
|Tombesi et al.| ( [2010a^ on the basis of their XMM-Newlon sample. 

In this work, we build upon our previous results to assess the 
properties of the disk-wind itself. The remainder of this paper is 
laid out as follows: in ^we outline the equations through which 
assess the wind properties, in p.2| we use the inferred values to per¬ 
form a qualitative assessment of the wind properties, before going 
on to conduct a quantitative correlation analysis of how the various 
wind properties scale with the AGN bolometric luminosity in p.3| 
In we then compare our results to those already present in the 
literature, discuss the wind acceleration mechanism which may be 
responsible for the observed outflows, and finally assess whether 
the observed wind is likely to be energetically significant in terms 
of feedback. We summarise our results and present our overall con¬ 
clusions in 


2 WIND PARAMETERS 

2.1 Projected distance 


The maximum column density Ah of gas along the line-of-sight 
(LOS) is given by Ah = n(r) dr, where n{r) is the average 
gas number density and r^ax is the maximum observed distance of 
the absorber. This, comb ined with the definition of the io nisation 
parameter ^ i Tarter, Tucker & Salpeter 1969| , yields 

the maximum distance that the absorber can be located from the 
ionising source given its observed column density and ionisation 
state. 


<^/Vh’ 


( 1 ) 


where Ljon is the source ionising luminosity integrated between 
1 — 1000 Rydberg. Conversely, a lower limit on r can be inferred by 
considering the escape radius of the gas given its observed velocity. 
For a simple Keplerian disc orbiting a black hole the escape veloc¬ 
ity at distance r is Vesc = In the limit that v* = Vesc, 

i.e., assuming that the measured outflow velocity along the LOS is 
equal to the escape velocity at observed radius r, we can set a lower 
limit on the location of the wind. 


^min 


2GMsh 


( 2 ) 


These relations hence allow upper/lower limits to be placed on the 
wind location given the measured parameters of the absorbing gas, 
albeit with large uncertainties. 


* Accessible at: http://heasarc.gsfc.nasa.gov/docs/archive.html 


2.2 Mass outflow rate 

The mass outflow rate of the wind, dM/dt = Afw, is a crucial pa¬ 
rameter and is the main means through which the overall flow en¬ 
ergetics are assessed. The mass outflow rate for an arbitrary wind is 
given by M* = A(r)p(r)v(r), where p(r) and v(r) are the density 
and velocity profile of the wind and A(r) is a factor which accounts 
for the geometry of the system. Assuming that the flow has con¬ 
stant terminal velocity v{r) = Vw and that the absorbing gas has 
cosmic elemental abundances (i.e., ~ 75% of its mass by hydrogen 
and ~ 25% by helium), p(r) ~ L2mpn(r) where nip is the proton 
mass and n{r) is the electron number density of the plasma. For a 
thin spherically symmetric isotropic wind A(r) = Qbr^, where the 
product Q.b 5/ 1 is known as the global filling factor and accounts 
for both the solid angle occupied by the flow (Q.) and how much 
of the flow volume is filled by gas (b). Thus, ~ Q.br^m^n{r)vvi 
where we have neglected the constant factor of order unity. 

The value b is extremely difficult to determine because 
it depends on the ionisation and dumpiness of the gas. At 
low-intermediate ionisation states the flow is likely to be 
clumpy/filamentary, while at high ionisation states it can be consid¬ 
ered largely smooth and of low density because the vast majority 

of elements are stripped of electrons. In the clumpy case, the col¬ 

umn density of the wind can be given by Ah ~ bn{r)8r, with b < 1 
implicitly allowing for inhomogeneities in the flow. Alternatively, 
at high ionisation states b ~ 1 and r —> r^ax- Substituting for rmin 
and Tmax then leads to algebraic upper and lower limits on the mass 
outflow rate: 

~ f2mpLion<^^'vw, (3a) 

~ 2flGMBHmpA h v„'. (3b) 

The major remaining uncertainty is the opening angle of the sys¬ 
tem Q.. Here, we adopt the average opening angle 12 = L6;r as in¬ 
ferred from the global detection fraction of Fe K winds reported in 
the literature (/ ~ 40%; |Tombesi et al.|2010a[[2014[ Paper 1). The 
assumption of a uniform wind geometry is clearly an oversimplifi¬ 
cation of real systems; in reality, numerous factors such as source 
luminosity, gas density, etc., will contribute to the shaping of an 
X-ray disk-wind and therefore the wind geometry is going to dif¬ 
fer on an object-by-object basis. Even so, we note that the average 
opening angle of £2 = L6;r is comparable to the wind opening mea¬ 
sured from the P-Cygni-like Fe K profile in PDS 456 (£2 is resolved 
to be > 2;r; see |Nardini et al.|2015^ . This suggests that it is a good 
approximation to the geometry of real disk-winds. 


2.3 Kinetic power 

Provided that the wind has already reached a steady terminal ve¬ 
locity by the point at which it is observed, the mechanical power 
imparted by expelling mass at a rate with velocity v„ is simply 
equal to its kinetic energy: = M^v^/2. We determine the range 

of likely kinetic power by substituting for Mw = 
leading to 

~ f2rapLion ^(4a) 

L™" ~ G£2mpMBHAHVw, (4b) 

Similarly, the rate at which the outflow transports momentum into 
the environment of the host galaxy is given by dpfdt = pout = 
■ Substituting again for leads to a plausible range between 


(c) ? RAS, MNRAS 000, ??-?? 
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Table 1. Summary of AGN and wind properties 


Source 

log Mbh 

log Lion 

log Ah 

log'? 

Vw/C 

(1) 

(2) 

(3) 

(4) 

(5) 


8.1 ±0.5 

44.7 

^■^•'-’-0.8 

4 45+0.40 
^•^^-0.46 

0 079+0 041 

U.U/Z_QQ3g 

3C 390.3 *'5 

Q q+0.2 

44.7 

>23.7 

>5.46 

0.145 ±0.007 

4C-t74.26‘® 

9.6*" 

47.0 

>21.8 

4.06 ±0.45 

0.185 ±0.026 

APM 08279 

lO.O* 

45.4 

23.0 ±0.1 

3 51+0.33 

0.285t!>:|i 

CBS 126*^ 

7.8T0.1 

44.2 

>23.7 

4 77+0-26 
' '-Q.il 

0.012 ±0.006 

ES0103-G035‘^’* 

7.4T0.1 

43.9 

>21.9 

4.36± 1.19 

0.056 ±0.025 

MCG-6-30-15'' 

6.7T0.2 

43.2 

22.2 ±0.1 

3.64±0.06 

0.007 ±0.002 

MR 2251-178® 

8.7T0.1 

45.2 

21.5±0.2 

3.26±0.12 

0.137 ±0.008 

Mrk279® 

7.5T0.1 

43.4 

23.4±0.3 

4 49+0.15 
^•^^-0.27 

0.220 ±0.006 

MrkVee** 


43.2 

99 7+0.2 
^^•'-0.3 

3.86l°:g 

U.U9y_QQ25 

NGC 1365*3.* 

7.6T0.7 

42.8 

22 7+^-2 

3.88±0.07 

<0.014 

NGC 3227* 

7.6T0.2 

42.5 

99 7+0.2 
^^•'-0.3 

'X S8+0-A2 

< 0.008 

NGC3516* 

7.6T0.1 

43.6 

22.6±0.2 

3.84±0.11 

0.004 ±0.002 

NGC 3783 * 

7.5T0.1 

43.6 

21.8±0.2 

3 48+0-*5 
■2-^°-0.07 

< 0.008 

NGC 4051* 

6.3 ±0.2 

42.5 

22.8 ±0.1 

■^•‘^^-0.50 

0.018+JJ;™ 

NGC 4151* 

7.1±0.1 

42.9 

>21.7 

3.69 ±0.64 

0.055 ±0.023 

NGC 4395 3 

4.7±0.2 

40.7 

22.8 ±0.3 

3.92±0.16 

< 0.001 

NGC5506®'’.*'’.* 

7.3±0.7 

43.7 

23.2±0.3 

^•^^-0.17 

0.246 ±0.006 

PDS456**.* 

9.4±0.3 

45.3 

23.0±0.1 

^•^°-0.15 

0.273 ±0.028 

SW 12127*3.* 

7.2+ 

43.7 

22.8 ±0.3 

4 16+®'^^ 

0.231 ±0.006 


Columns:- (1) SMBH mass, in units of M©; (2) absorption-corrected 1-1000 Rydberg 
ionising luminosity (in units of ergs“^), as determined from the best-fit continuum model 
outlined in Paper I; (3) measured wind column density, in units of cm“^; (4) wind ionisation 
parameter, in units of erg cm s“^; (5) inferred wind outflow velocity, as measured in Paper I. 


References:- jPeterson et al.| |2UU4|; jUdn et al.| |2U12|; jChatteriee et al.M2Ull>: 
Tombesi et al.H2012b); ^Bentz et al.H20Q9a); ^Khorunzhev et al.H2012|i; ^McHardy et al 


^Peterson et al. ] j2004); ^Edri et al.| j2012); ^Chatterjee et al.|j201lt; 


I 2005t; ^Bentz et al.li 2Q09bt; iPapadakisH2004|; ^Nikotajuk Czerny & 
^ Reeves et al. j2Q09 ; ^ jMalizia et al H2Q08t ; ^^isaliti et al.H2007t; 


Czerny et al 

Grupe et ^ 2004f ; ^ ^Woo & UrryH200'^ ‘ jSaez, Chartas & Brandt| 2009) . 


Gu^n 


2009 : 


2001 


Notes:- *Mbh is taken as the mean of extreme values found in the literature; 'En'ors 
on Mbh not present in literature; ^Full designation APM 08279+5255; ^Full designation 
SWIFT J2127.4+5654. Where a source has an FeK wind detected in more than one obser¬ 
vation, the mean Ah, ^ and Vw values are reported here; see Paper I for further details. 


~ fimpLion^ (5a) 

~ GflmpMfiHVH, (5b) 

2.4 Other parameters 

2.4.1 Black hole masses 

Our method of estimating the inner radius of the outflow, is 
proportional to the escape radius of the black hole, resc, which in 
turn is proportional to the mass of the central black hole (BH), Afgu- 
By extension, this also means that the other lower limiting quanti¬ 
ties in equations l |3b^ , l |4b[ l and l |5b^ are also proportional to 
Mbh- Estimates for gathered from the literature are collated 
in Table [T] We obtain most of the Mbh values through the numer¬ 
ous reverberation mapping studies available in the literature (i.e., 
[Peterson et al.|2004l [Bentz et al.|2009a|b| [Edri et al.|2012^ which 
tend to offer relatively tight constraints on the mass of the central 
object; in the cases where an AGN had been subject to reverbera¬ 
tion mapping multiple times, we adopt the most recent estimate on 


Mbh only. There are several empirically determined mass estimates 
for the BHs in 3C 111, ESO 103-G035, NGC 1365, NGC 5506 and 
PDS 456. For these AGN we report the mean value, with the asso¬ 
ciated errors taken as half of the range between the minimum and 
maximum values to account for the uncertainty in the individual 
estimates on Mbh- For CBS 126, where we were unable to find any 
robust mass estimates in the literature, we infer Mbh from the BLR 
line width and luminosity scaling relation ofJKaspi^^^^lJ d2000|l 
using the appropriate spectral values listed in jGtupe et al.|l ^04|>. 
For 3C 390.3 we use the range of values inferred hy Dietrich et al.j 
\2Q\2) . In all other sources we take Mbh directly from the refer¬ 
ence listed in the table footnote. The sample encompasses almost 6 
orders of magnitude in black hole mass, with logMnn ranging from 
~ 4.7 in the dwarf-Seyfert galaxy NGC 4395 jEdri et al.|2012| l, all 
the way up to an estimated ~ 10 in the high luminosity BAL quasar 
APM 08279-1-5255 jSaez, Chartas & Brandt|200^ . 


© ? RAS, MNRAS 000, ??-?? 
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Table 2. Summary of inferred wind parameters. 


Source 

logrmm 

(cm) 

logrmax 

(cm) 

logM™" 

(gs^‘) 

logM™ 

(gs^') 

logL™" 

(ergs^*) 

logL““ 

(ergs^') 

3C111 

15.9±0.1 

17.3±0.6 

24.9 ±0.1 

26.4 ±0.6 

43.3 ±0.1 

44.7 ±0.6 

3C 390.3 

15.7±0.1 

< 16.7 

>25.5 

<26.6 

>44.5 

<45.6 

4C-1-74.26 

16.5±0.1 

< 20.0 

>24.8 

27.9±0.5 

>44.0 

47 1+0-5 

APM 08279 

16.6±0.1 

19.8±0.3 

26.3 ±0.1 

29.6 ±0.3 

45.9 ±0.1 

49.1 ±0.3 

CBS 126 

- 

< 15.9 

- 

24.8±0.3 

- 

41.6±0.2 

ESO 103-G035 

15.4+« 

< 18.8 

>23.0 

25.6±1.2 

>41.1 

43.7 ±1.2 

MCG-6-30-15 

16.5 ±0.3 

17.4±0.1 

23.8±0.3 

24.7 ±0.1 

40.1 ±0.3 

41.0±0.1 

MR 2251-178 

15.9±0.1 

20.4 ±0.3 

24.3 ±0.1 

27.8±0.7 

43.3 ±0.1 

46.7 ±0.7 

Mrk279 

14.3±0.1 

15.6±0.5 

24.3 ±0.1 

25.6±0.5 

43.6 ±0.1 

44.9 ±0.5 

Mrk766 

14.5 ±0.1 

16.7±0.3 

23.1 ±0.1 

25.2±0.3 

40.9 ±0.1 

43.1 ±0.3 

NGC 1365 

- 

15.2±0.1 

- 

24.1 ±0.1 

- 

40.5 ±0.1 

NGC 3227 

- 

16.0±0.4 

- 

23.6±0.4 

- 

39.7 ±0.4 

NGC 3516 

- 

17.2±0.3 

- 

24.6 ±0.3 

- 

40.5 ±0.3 

NGC 3783 

> 17.3 

18.4±0.3 

>24.2 

<25.5 

<40.5 

<41.9 

NGC 4051 

15.3±0.1 

15.7±0.2 

23.6±0.1 

24.0 ±0.2 

40.7 ±0.1 

41.2±0.2 

NGC4I51 

15.1 ±0.4 

< 18.1 

>22.6 

25.2±0.6 

>40.7 

43.4±0.6 

NGC 4395 

- 

13.9±0.4 

- 

<21.5 

- 

<36.1 

NGC 5506 

14.0 ±0.1 

15.4±0.5 

23.9±0.1 

25.3 ±0.5 

43.3 ±0.1 

44.7 ±0.5 

PDS 456 

16.0±0.1 

18.2±0.2 

25.8±0.1 

27.9 ±0.2 

45.3 ±0.1 

47.5 ±0.2 

SWJ2127 

13.9±0.1 

16.8±0.5 

23.4±0.1 

26.2 ±0.5 

42.8 ±0.1 

45.6±0.5 


2.4.2 Bolometric luminosity and Eddington ratio 


3 DATA ANALYSIS 


We estimate the bolometric luminosity for each AGN: Lbol = 
I‘’boli2-l0k;ev> where Li-lOkeV is the unattenuated source luminos¬ 
ity integrated between 2— lOkeV and iCbol is the bolometric cor¬ 
rection factor. Various studies have shown that the spread of 
amongst individual AGN is quite large (e.g., [Elvis et al.|[T9M| ), 
whilst the correction appropriate for a particular AGN can be a 
function of luminosity ([Marconi & H unt|2003[ [Hopkins, Richard sj 
& Hemquist 2007 Eddington ratio i Vasudevan & Fabian[[2007[ 
Lusso et al.|2010[l, or both jLusso et^ 2010^, such that the uniform 


application of a single bolometric correction factor may be inappro¬ 
priate for a heterogeneously selected sample. We therefore primar¬ 
ily use the values listed in the works of Vasudevan & Fabiar^ 
( [2007t ; [Vasudevan et al.|P00^ and [Vasudevan et al.[ ( [2010) which 
are empirically determined on the basis of the broad-band spectral 
energy distributions of each AGN. Several of the AGN in the sam¬ 
ple are not listed in these works. For these sources, we searched 
the literature for an appropriate bolometric luminosity and used 
that instead (see Table [^caption). Only for 3C 111 and CBS 126 
were we unable to locate either a bolometric correction factor or 
an definite and empirically measured bolometric luminosity; for 
these two sources we simply assume that )q,oi ~ 30, which is simi¬ 
lar to that found previously for other sources which harbour UFOs 
(e.g., [Tombesi et al.|2Q12a[ Values forL2_iokeV (hereafter denoted 
Lx), Kbol and the resultant estimate of Lboi are all listed in Ta- 
ble[^ For completeness, we also compute the likely Eddington ratio 
of each AGN: A = Lboi/Leddi where L^dd = ^KGmpM^-nCdj^ ~ 
1.26X 10^*(Mbh/Mq). 


3.1 Preparing the data 

Several of the AGN in the have Fe K absorption detected in multi¬ 
ple epochs (see Paper I), such as is the case for 3C 111, Mrk766, 
NGC 1365, NGC 3227, NGC 3783, NGC 4051, or have more than 
one absorption trough which comprises a multi-velocity system 
(PDS456 and APM 08279-1-5255). For these AGN considering the 
mean parameter value and accounting for intrinsic outflow variabil¬ 
ity between epochs by folding the uncertainties on each of the indi¬ 
vidual measurements into the error bar. In the majority of cases we 
find the individual measurements to be largely consistent within the 
errors, suggesting that the outflow is persistent across the different 
epochs. Only in 3C 111, Mrk766 and NGC4051 do there appear 
to be significant differences at the 90% level (see Table[TJ, and this 
is reflected in the broad error bars for these sources. For PDS 456 
and APM 08279-1-5255, which both appear to have two FeK ab¬ 
sorption systems at different velocities, we adopt the mean velocity 
and again fold the range of possible velocities into the error bar. 
Table[T]summarises the measured parameters of the winds detected 
in the Suzaku sample. 


3.2 Parameter constraints 

Using the values in Table we first computed the wind parame¬ 
ters in standard units and then normalised them to the appropriate 
values for a given black hole mass; the resultant values are noted 
in Tables and[^ respectively.We normalised the specific parame¬ 
ters as follows: (i) the distance (rmax, rmin) to units of Schwarzchild 
radius (ii) the kinetic luminosity (L^^jL™") to 

A-Edd. (iii) the mass outflow rate to the Eddington 

accretion rate A^Edd = ^Edd/tic^ (assuming t] = 0.06 for the ac¬ 
cretion efficiency), and (iv) the momentum rate to the Eddington 


(c) ? RAS, MNRAS 000, ??-?? 
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momentum rate />Edd = ^Edd/c- Algebraically, this yields a set of 8 


normalised equations: 

^max/^s ~ i (6a) 

^mm/''s = (c/vw)^, (6b) 

^w“/^Edd = (4;rG<^MBHne)^‘fiCTTnHAonVwI]C, (6c) 

= (2;rvwne)-lflcTT«HAfHJ]c, (6d) 

L™“/LEdd = (47tGMBHCi^ne) 'fiCTxUHGionVw: (6e) 

= (4;rc«e)^'fiO’T«H%Vw, (6f) 

/>w“/PEdd = (47rGMBH<^ne)^'f20’T«HAonVw, (6g) 

Pr7PEdd = (4;rne)^‘f2cTTAH. (6h) 


Errors on the normalised parameters were determined by both prop¬ 
agating the uncertainties on the measured XSTAR values through 
the various equations, and also taking into account the error on Mbh 
where applicable. For CBS 126, NGC 1365, NGC 3227, NGC 3516 
and NGC 4395, which all have relatively slow wind velocities of the 
order Vw < 0.01 c (see Table^, we found that > rmax (within 
the errors) which indicates that either: (i) the outflow may not have 
achieved the requisite escape speed for its observed location, in 
which case the wind material may not escape into the host galaxy 
unless it is subject to an additional acceleration mechanism once it 
has crossed the LOS, or (ii) that the wind steam line flows almost 
perpendicularly along the LOS such that we are only seeing the tan¬ 
gential component of v*. While the associated error bars in these 
slow sources may overlap (suggesting that the upper and lower lim¬ 
its may be formally consistent) we choose to only report the upper 
limiting quantities for these five AGN in Tables and so that 
we do not introduce an artificially tight (and potentially mislead¬ 
ing) constraint on the properties of the wind. In the remainder of 
our analysis we separate the outflow sample into two groups based 
on their velocity: those with Vw ^ 0.01 c (3000km s^*) are here¬ 
after classified as ‘slow’, whilst those with Vw > 0.01 c are ‘fast’. 
Whilst the chosen velocity threshold is essentially arbitrary, this 
classification scheme provides a useful means to distinguish be¬ 
tween those pristine ‘disc-winds’ from those which may just be the 
higher-ionisation component of the inner-BLR or the more distant 
warm absorber, with which they probably share an overlap in ve¬ 
locity space and location (e.g., see |Tombesi et al.|2013| l. 

In Figure[2we show constraint diagrams for the various wind 
parameter^ From Figure [^-b it is clear that the absorbing ma¬ 
terial is distributed across several orders of magnitude in distance 
from their black hole, with r spanning ~ 10*^^^® cm or ~ 10^^® Xj. 
By considering the mean rang^we see that the fast outflows ap¬ 
pear to be located between (r) ~ 10*^^^^ cm (~ 10^^“^ rs) dis¬ 
tance range. This corresponds to ~ 0.0003 — 0.3 pc and implies that 
whilst there may be some limited overlap between the fast FeK 
outflows and the traditional soft X-ray warm absorber (which is 
usually inferred to be on parsec-scale distances, e.g., Blustin et al.| 


|2005[[Kaastra et al.|2012[[Crenshaw, Kraemer & George|2003 and 

references therein), the former are generally located much closer 


^ Note that we do not show a constraint diagram for because it is simply 
proportional \aP,/v^ and therefore follows a similar overall distribution to 
the one for L„ 

^ Here, ‘mean range’ refers to the range spanning between the means for 
the upper and lower limit values, respectively, i.e., it corresponds to the 
range between between (rn,ax) —> (rmin) for the fast outflows (see caption 
for more details) 


to their central black hole. This is in agreement with other re¬ 
sults in the literature which find the highly-ionised Fe K outflows 
most likely originate in a wind from the inner regions of the ac- 


cretion disc (e.g., Proga & Kallman|2004||Schurch, Done & Proga| 

20091[Sim et al.|2008 

King 2010l, and is also consistent with what 

Tombesi et al.[([2012a 

found on the basis of the XMM-Newton out- 


flow sample. 

Following our estimates on r. Figure [T]; shows an anal¬ 
ogous plot for Mw The constraints on also vary signifi¬ 
cantly, spanning almost 8 orders of magnitude, and ranging from 
< 10^^ gs^* in the dwarf Seyfert galaxy NGC 4395 all the way 
up to ~ lO^'^gs^* in the massive quasar APM 08279-1-5255. The 
other winds have estimates distributed between these two ex¬ 
tremes, with the mean range for the fast outflows falling between 
(Afw) ~ iq24-26 gg-l Q Qj _ IMq yr). Interestingly, and while 
there remains considerable uncertainty in each case, (M™“) is of 
the order of the Eddington rate, i.e., ~ MEdd> while even 

(M™”) is still ~ 1% of Medd- This immediately suggests that the 
winds are transporting a substantial amount of material into their 
host galaxies. The ejected mass has a large kinetic power: Fig¬ 
ure shows that while the range of L„ again spans several orders 
of magnitude, the overall mean for the fast flows is constrained 
to between (Lw) ~ 10^^ and 10^^ erg s^^. The normalised val¬ 
ues show that, on average, (L„) is < LEdd> with a mean range of 
(Lw) ~ (0.1 — 10%)LE(id- The wind mass outflow rate and the as¬ 
sociated kinetic power are therefore significant fractions of the Ed¬ 
dington limited values. In §4.3| we extend these results to assess 
whether the wind is energetically significant in terms of feedback. 

As mentioned previously, our estimates on r are subject to 
large uncertainties. It can therefore be useful to compare a few of 
our inferred values with some more robust measurements available 
in the literature. We consider here the exemplar cases of PDS 456 
and Mrk766. Currently, PDS 456 is the only AGN where the lo¬ 
cation of the Fe K-shell wind has been robustly constrained from 
discrete line variability. On the basis of a recent Suzaku campaign 
[Gofford et al.| ( |2014| l constrained the Fe K absorber in PDS 456 
to ~ 100— 1700^5. This is entirely consistent with what we in¬ 
fer using equations jnuL Similar is also true for Mrk766:| Miller| 
|et al.|p()07^ and |TumM et al.| ( |2007^ estimates the absorbing gas 
to be located on distances of the order 30 — 50 Xj, while |^saliti| 
|et al.| ( [^ll[ l later showed that both the broad-band spectral vari¬ 
ability and the Fe K absorption could also be attributed to a strat¬ 
ified, inhomogenous and clumpy absorber which occults the AGN 
at a distance of ~ 500 — 5000 r^. These estimates are again consis¬ 
tent with what we find here (see Figure [TJi-b), suggesting that the 
two studies are perhaps probing the same layer of gas. Elsewhere 
in the literature, [Miller et al.| pOlOj l used a detailed spectral tim¬ 
ing analysis to de-convolve the X-ray spectrum of NGC 4051 into 
its constituent absorption- and reflection-dominated components, 
with the hard X-ray reverberation signal suggesting that the reflect¬ 
ing gas is ~ 100 — 600 from the black hole (at the 90% level, see 
[Miller et al.|2010] l. This is around one order of magnitude closer 
than where we estimate the Fe K absorbing gas to be in this source. 
One possibility is that the reflecting layer detected by [Miller et al.[ 
|2010[is physically distinct to the FeK absorber, possibility asso¬ 
ciated with the Compton-thick occluding clouds posited by [Tatuin| 
[et al.[ ( [20T3] 2014, in prep). 


3.3 Correlations with bolometric luminosity 

Supplementing the constraint diagrams we also carried out a cor¬ 
relation analysis. This analysis bears notable similarities to the one 
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Table 3. Summary of normalised wind parameters. 


Source 

^bol 

^bol 

logr^in 

('•.) 

logfmax 

('•.) 

logM™" 

{Mem) 

logM“‘*’‘ 

(Vedd) 


log!™" 

(f-Edd) 

logL““ 

(f-Edd) 

3C111 

44.35±0.18 

45.7 

20 

2.3±0.1 

3.7 ±0.6 

-1.3±0.1 

0.1±0 

6 

-2.9±0.1 

-1.5 

±0.6 

3C 390.3 

44.42 ±0.01 

45.5 

13.3" 

I.7±0.1 

<2.8 

> -1.1 

< -0 

0 

> -2.1 

< 

-2.8 

4C+74.26 

44.95 ±0.01 

46.3 

? 

1.5±0.1 

<5.0 

> -2.8 

0.2 ±0 

5 

> -3.8 

-0.6 

±0.5 

APM 08279 

46.73 ±0.03 

47.4 

_1 

I.l±0.1 

4.4±0.3 

-1.9±0.1 

1.4±0 

3 

-2.3±0.1 

1.0±0.3 

CBS 126 

43.82±0.02 

45.1 

20’ 

- 

<2.6 

- 

-1 1+® 

i-i_0 

2 

3 

- 

-4 

1 + 

p p 

ESO 103-G035 

43.47 ±0.01 

44.6 

12.6‘' 


<6.0 

>2.6 

0.0 ±1 

1 

> -4.4 

-1.8 

±1.1 

MCG-6-30-15 

42.77 ±0.01 

44.0 

16.1" 

4 3+^-^ 
^•■^-0.2 

5.23±0.1 

_1 1+0.3 
^•^-0.2 

-0.2±0 

1 

1 

1 + 

p p 

-3.8 

±0.1 

MR 2251-178 

44.61 ±0.01 

46.0 

22.0" 

1.7±0.1 

5.2±0.1 

-2.53±0.1 

1.01° 

5 

8 

-3.55±0.1 

0 

1+0.5 

•-^-0.8 

Mrk279 

42.78 ±0.01 

44.0 

10. l" 

1.3±0.1 

2 6+^'^ 
"^•^-0,4 

-1.36±0.1 

-0.11° 

6 

4 

-1.98±0.1 

-0 

n+Q.6 

•'-0.4 

Mrk766 

42.73 ±0.05 

44.1 

48.9" 

2.8±0.1 

5.0±0.3 

-1.3±0.1 

0.9 ±0 

3 

-3.4±0.1 

-1.2 

±0.3 

NGC 1365 

42.32 ±0.02 

44.1 

101.6" 

- 

2.1±0.1 

- 

-1.7±0 

1 

- 

-5.2 

±0.1 

NGC 3227 

42.13 ±0.06 

43.3 

15.3" 

- 

2.9±0.3 

—0 4+^-'* 

-2.2±0 

3 

- 

-6.1 

±0.3 

NGC 3516 

43.07 ±0.01 

44.1 

14.0" 

- 

4.1±0.3 

- 

-1.13±0 

3 

- 

-5.23 

±0.3 

NGC 3783 

43.14±0.05 

44.5 

23.7" 

>4.3 

5.4±0.3 

> -1.5 

< -0 

1 

<-5.1 

<-4.3 

NGC 4051 

41.65 ±0.04 

43.2 

41.1‘' 

3.5±0.1 

3.9±0.2 

-0.9 ±0.1 

-0.4 ±0 

2 

-3.7±0.1 

-3.2 

±0.2 

NGC4151 

42.34 ±0.02 

43.6 

18.3‘' 

2 5“^^'^ 
^■■^-0.3 

<5.6 

>-2.7 

0.0 ±0 

6 

> -4.5 

-1.8 

±0.6 

NGC 4395 

40.38 ±0.01 

41.7 

22.5" 

- 

3.8 ±0.4 

- 

< -1 

4 

- 

< -7.5 

NGC 5506 

43.19±0.01 

44.4 

16.8" 

I.2±0.1 

2.6±0.5 

-1.6±0.1 

-0.2±0 

5 

-2.1 ±0.1 

-0.7 

±0.5 

PDS456 

44.62 ±0.01 

47.0 

_2 

1.1±0.1 

3.3±0.2 

-1.8±0.1 

0.4 ±0 

2 

-2.2±0.1 

0.0 ±0.2 

SWJ2127 

43.15±0.01 

44.5 

7 

I.3±0.1 

4.1±0.5 

-2.0±0.1 

0.8±0 

5 

-2.6±0.1 

0.3 

±0.5 

Notes:- 'The bo 

lometric luminosity of APM 08279+5255 is take 
uming a conservative magnification factor of fir 

n as Lhoi = 7 X 10'"u, ' Lp, 

Lewis et al. 

1998 

Riechers et al 

2009 

Saez & 

|Chartas|201 1|, as^ 

= 100 ^Egami 

et al.|2000k; for PDS456 

taken from Reeves & 

: Turner 


the r — Ifboi coiTelation of |Zhou & ZhaoH2010^ using the best-fit F values that we found in Paper I. ^Lx here corresponds to the unattenuated 
source luminosity between 2 — lOkeV. 


recently conducted by |Tombesi et al.| ( |2013^ for the XMM-Newton 
outflow sample. There, the authors combined the parameters of the 
hard-band absorber, as measured in their FeK absorption survey 
with XMM-Newton, with those for the soft X-ray warm absorber 
collated from the literature to determine how the AGN wind pa¬ 
rameters vary globally in relation to their distance from the black 
hole ( [Tombesi et al.|20T^ . Here, noting that Suzaku does not have 
the soft X-ray energy resolution necessary to constrain the veloc¬ 
ity of the soft X-ray absorbing gas, we consider only the Fe K ab¬ 
sorbers and do not consider the warm absorber. Whilst this leads 
to a relatively small dynamic range in terms of their distance from 
the black hole, the heterogeneous nature of the Suzaku sample en¬ 
compasses a significantly broader range of AGN in terms of their 
bolometric luminosity (Tbol)- The Suzaku sample therefore offers a 
useful opportunity to probe how the wind parameters may vary in 
comparison to the luminosity of their host AGN. 

We compute linear regressions of the form log(y) = a + 
j8 log(x), where a and j3 are the intercept and slope of the straight 
line fit, re spectively, u sing the lirmiix_err Bayesian regression al- 
gorithrr|^ i Kelly|2Q07 i. This routine employs Markov Chain Monte 
Carlo (MCMC) techniques to self-consistently account for mea¬ 
surement errors and intrinsic scatter in the data, whilst also al¬ 
lowing for limited censorship (upper-limits only) in the indepen¬ 
dent variable. We computed regression parameters from 100,000 
MCMC realisations of the data by the linmix_err routine. This 
yields posterior distributions for the intercept a, slope j8, standard 


available in IDL from: 


deviation (scatter) of the data a, and the Pearson linear correlation 
coefficient Rp. We adopt the median value of the posterior distribu¬ 
tion as our ‘best-fit’ to the data, and estimate the parameter errors 
as the range of simulated values which encompass 90% of the pos¬ 
terior distribution about the median. 

We note the slow winds have Vw ^ 3000 km s^ * which is sim¬ 
ilar to that measured in the traditional soft X-ray warm absorber. 
It is therefore unclear whether these winds are associated with a 
bona fide disk-wind or alternatively with another layer of gas which 
is more distant from the black hole, e.g., with a higher ionisation 
component of the warm absorber. For this reason, we fit all subse¬ 
quent regression analyses to the fast winds only, Vw > 3000 km s^^. 
Owing to their high velocities, these faster winds are kinematically 
distinct from the soft X-ray absorbing gas and thus may more ob¬ 
viously represent the signature of a pristine disk-winds. This limits 
our analysis to the winds which are capable of escaping the gravi¬ 
tational potential of the central black hole (based on our LOS), and 
are therefore more likely to have an effect on the host galaxy in 
terms of feedback. 

In Figure|^we show how log(AH/cm^^), log(i^/ergcms^*) 
and log(vw/kms^*), as measured in Paper I and summarised in 
Tablevary with Lboi- From Figures]^ and[^ we find no dis¬ 
cernible relationships between either Ah or ^ and Lboi (/3 = 0 at 
90% confidence): the log(AH/cm^^) values cover a wide range for 
a given value of Lbob while log(^/ergcms^*) appears to cluster 
at around ~ 4, mirroring the median of the parameter distributions 
that we found in Paper I. This also appears to be true for both the 
fast and the slow wind sub-samples. Conversely, there is a correla¬ 
tion present between Vw and Lbob with the faster winds tending to 
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Figure 1. Constraint diagram showing the parameter ranges occupied by the 5Mz^?/:M-detected Fe K outflows. The left hand side panels shows the raw limiting 
values (taken from Table[^ and the right hand side shows the normalised ones (from Table[^. The top, middle and bottom rows are the constraints on distance, 
mass outflow rate and kinetic power. In all panels, the red and blue lines show sources with > 3000kms“^ and Vw < 3000km s“^, respectively. For Mrk766, 
NGC4051 and PDS456, the grey boxes show the range in location measured by |Risaliti et al.|(201 l^ ; [Miller et al^j2007) ; |Tumer et al.H2007) ; |Miller et al.| 
j2010) and |Gofford et al.]j2014) on the basis of variability or X-ray reverberation studies, respectively (see text for further details). The blue/shaded areas 
show the ‘mean range’, spanning the mean of all upper and all lower limiting values of the fast winds only. 


be being observed in more luminous AGN. The fast winds have a 
constrained slope of p = 0.4^q 2, with Vw and jS 0.5. In 

contrast, the slower systems, i.e., those with < 0.01 c, appear to 
be isolated to the lower left quarter of the plot which suggests that 
low velocity winds are preferentially located in lower luminosity 
systems. This is consistent with the overall picture that wind veloc¬ 
ity is correlated with L^oi; fitting a correlation to the entire sample 
yields j3aii = 0.5 ^q 2. This is formally consistent within the errors 
but does suggest that the inclusion of the slow winds may skew the 
relationship to a slightly steeper slope. 

Figure then shows how the raw wind parameters vary with 
Lboj. Here, and in subsequent figures, each data point corresponds 
to the mid-point in the range for each parameter, calculated as, e.g.,: 

_ (rmax + perr) + (rn,in-nerr) 

^mid 2 ’ kt/ 

where nerr and perr are the negative and positive error, respec¬ 
tively. The associated error bars then denote the range between the 
maximum and minimum values, including the uncertainties aris¬ 
ing from the input variables. In equation|^we are essentially only 


considering the upper- and lower-limits for all derived values in 
the sample, such that in the cases where the minimum and maxi¬ 
mum values are themselves unconstrained (i.e., only and upper- or 
lower-limit is available with no associated error) we are able to use 
the censored value when calculating the mid-point. This ensures 
that all of the available data is included in our analysis. 

To first order, we find that all of the raw wind parameters are 
plausibly correlated with Lbol to some degree. Whilst the slope for 
Lbol “ f is tentative, and marginally consistent with j3 = 0 at the 
90% level (Figure]^ the overall distribution is skewed towards a 
positive relationship among the data which suggests that the winds 
are observed at larger radii in brighter AGN, as would be expected 
given the larger typical size-scales in these systems. Figure 
then shows that there is a positive relationship between Lbol ^ 

(j3 = O.O^og), such that the winds in more luminous AGN are cor¬ 
respondingly more massive. In both of these cases the slow winds 
appear to follow the same overall trend as the fast ones, despite not 
formally being included in the regression computation. Similar is 
also true for L(,oi — and Lboi — Pw. which share strong positive 
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Figure 2. Scatter plots showing logLboi versus: (a) log(A^H/cni“^), (b) 
log(^/ergcms“^) and (c) log(vw/knis~*), as noted in Tabled The red 
and blue data-points correspond to AGN with Vw > 3000 km s“~and < 
3000 km s“ ^, respectively. The solid black line corresponds to the ‘best-fit’ 
linear regression to the fast winds, as estimated from 100,000 MCMC reali¬ 
sations of the dataset using the liniiiix_err Bayesian regression algorithm 
jKelly|2007) , while the grey shaded area denotes the region which contains 
90% of the posterior probability distribution. The Parameters of the fits are 
reported in Table 

slopes of P = 1.5+1^;^ (R^ ^ 0.54) and j3 = 1.2+°;^ (Rp > 0.55), re¬ 
spectively. The observed slopes here likely stem from their mutual 
dependence on vv: Tw and pw v^. Indeed, the Lw 
and pw ^o'l ** relationships implied here are formally consistent 
with the v„ oc found earlier in Figurej^. This ultimately sug¬ 
gests that the wind velocity may be the driving factor in the ob¬ 
served relationships, with the uncertainties on the other wind pa¬ 
rameters only adding to the underlying scatter of the data. In con¬ 
clusion, these results indicate that disk-winds in higher luminos¬ 
ity sources are (tentatively) located farther away from their cen¬ 
tral black hole in terms of absolute distance (as expected for larger 
Mbh). they contain more material, and they are accordingly more 


Table 4. Summary of linear regression parameters: log(y) = a + /I log(x) 


Relation 

a 

P 

a 

Rp 

log(x) log(y) 

(1) 

(2) 

(3) 

(4) 

Tbol Al 

21-illU:? 

U.U_02 

0.6l!>;i 

0 11+039 
^•^^-0.65 


0.3^’:* 

^•^-0.2 

0.41“ 

^•■^^-0.60 

Vw 

-13 6+^2-^ 

0.41“ 

o-+o:i 

0 fSQ+0'26 

'-'•oy-0.47 

r 

— 8 4+38-3 

0.61“ 

0.81“ 

0 71+0.28 
^•'^-0.93 


—13 

0 9+0-8 

^•^-0.4 

>0.34 

Z/w 

973 c+23.6 
^•^•'^-44.7 


^•^-0.7 

>0.54 

Pw 

1 0 1 +28.8 

^°-^-36.9 

1 9+0.8 
^•^-0.7 


>0.55 

rju 

24 1 +^4-® 
^^•-^-34.3 

'-'•^-0.8 

o-slgl 

'-'•O'-'-0.38 

^w / ^Edd 

0 4+29-5 
^•^-29.6 

U.U_o7 

0-5l!;;I 

0 04+^'^® 

^•^^_0,99 

•Zw / ^Edd 

— 27 

^+^-35.2 

^•^-0.8 

o+ISI 

0.751“^ 

Pw / pEdd 

— 12 4+30-5 
^+•^-30.9 


^•^-0.4 

0.56l«:« 

Pbol Pw 

4+21.6 

■^•^-27.9 

1 9+0.8 
^•^-0.7 

0 7+^'^ 
^+-0.5 

>0.55 


Notes: All regressions were performed using the linmix.err Bayesian 
regression routine jKelly|2007) which takes into account both measurement 
and upper limits in the y-variable. The noted best-fit values correspond to 
the median value of the Posterior probability distribution, as simulated from 
100,000 MCMC realisations of the data, with errors taken as the range of 
values which encompass 90% of the Posterior probability. See |Kell^ j2007^ 
for more details. Nominal values correspond to a fit to the entire sample, 
while those in brackets are for an analogous fit with NGC4395 removed. 
“Regression computed with lower limits excluded. 

Columns: (1) slope of the linear regression; (2) normalisation/intercept of 
the best-fit regression line; (3) standard deviation of the intrinsic scatter in 
the data; (4) best-fit linear correlation coefficient. ±1 denote perfect posi¬ 
tive/negative correlations, respectively. 

energetic. The driving factor behind these relationships appears to 
be the wind velocity, which is faster in more luminous objects. 

Similar scatter plots for the normalised wind parameters are 
shown in Figure]^ In contrast to before, this time we are unable to 
discern any linear correlations in the normalised Lt,oi — (t/ts) and 
^bol ^ (^w/4^Edd) plots. In fact, and whilst there is a weak negative 
correlation suggested in the case of the former, both of the relation¬ 
ships are formally consistent with P =Rp = 0 such that the r/rs and 
Mw/tk?Edd ratios are largely independent of Lbol- The fast winds are 
therefore observed at similar relative distances from the SMBH re¬ 
gardless of the source luminosity, and their mass outflow rate is 
a similar fraction of MEdd- The slow winds also appear to follow 
consistent relationships. Both L„/LEdd (Figure]^) and Pw/pEdd 
(Figure]^) are also formally consistent with a slope of zero, al¬ 
though weak positive correlations are suggested. In both Figure]^ 
and S’ the slow winds appear to lie below the computed regres¬ 
sion line for the fast winds. Their location in these plots appears to 
mirror that seen in Figure]^ for the Lboi — Vw correlation, which is 
again consistent with the wind velocity being the parameter behind 
the relationships. 

For completeness, we also tried to determine whether the pa¬ 
rameters of the wind scale with the Eddington ratio of the AGN. 
Unfortunately, however, the relatively tight clustering in the data at 
X = 0.1, compounded by the lack of Fe K wind parameters in low 
Eddington ratio sources (i.e., X < 10^^) sources, meant that we 
were unable to discern any plausible relationships among the data. 

As a consistency check, we also checked for evidence of relation¬ 
ships amongst the raw Lx values finding that the same correlations 
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are still present in the data, albeit at lower significance level. The 
presence of the same correlations in both the contiguous Lx and the 
non-simultaneous Lbol indicates that the correlations are not driven 
by uncertainties in fCboj. 


4 DISCUSSION 

4.1 Comparisons with previous work 

4.1.1 The XMM-Newton outflow sample 

In Paper I we showed that ~ 40% of the AGN in the Suzaku 
sample exhibit evidence for highly-ionised winds in their FeK 
band. This is in agreement with the detection fraction reported by 
[Tombesi et al.| ( |2010at for XMM-Newton, and corroborates their 
conclusion that such winds are either (i) persistent over the ac¬ 
tive phase and thus have a large covering fraction, or (ii) tran¬ 
sient in nature and have a covering fraction of unity for only a 
fraction of the time. We also showed that the log(AlH/cm^^) and 
log(i^ /ergcms^^) distributions are also entirely consistent between 
the two samples, covering ranges of 21 < log(AH/cm -^) s: 24 
and 2 < log(<^/ergcms^*) ^ 6 in both, respectively, with the out¬ 
flow velocity Vw also being largely consistent on the basis of a 
Kolmogorov-Smimov (KS) test. By comparing the distributions 
of wind parameters found in this work to those determined by 
|Tombesi et al.|j2010a| > we similarly find that the FeK absorbers 
detected in both the Suzaku and XMM-Newton studies are typ¬ 
ified by the same range of physical parameters and occupy the 
same parameter space in terms of their overall location and energet¬ 
ics: 1 < log{r/r,) < 3 (~ 0.001 - 1 pc), -2 < log(Mw/MEdd) < 0 
(~ 0.001 — 10 Mq yr) and —3 < log(Lw/TEdd) iS 0- Their medians 
are also similar with log(r/rs) ~ 16(17), log(Mw/gs^*) ~ 25(25) 
and log(Lw/erg s^^) ~ 44(44) fov Suzaku (XMM-Newton), respec¬ 
tively, in natural units. 

The Suzaku outflow sample therefore robustly confirms that 
FeK absorption are a real physical component of emergent X- 
ray spectrum, that the implied winds likely have a large covering 
fraction, and that they are typified by large column densities and 
high ionisation. The range of outflow velocities, which are typi¬ 
cally Vw is 10,000 km s^ *, but can be as low as < 1,000 km s^ ^ in 
some cases, also suggests that the highly-ionised winds share an 
overlap in velocity space with the traditional soft X-ray warm ab¬ 
sorber. This is in line with the results of [Tombesi et al.||2013) who 
argue that some of the soft X-ray warm absorbers could be the arte¬ 
facts of accretion disc winds which have propagated farther away 
from the SMBH. 


4.1.2 Comparison with King et al. (2013) 

Another interesting study that we can consider is that of [King et al.j 
( [2013^ , who compare the kinetic power of both the warm absorber 
and relativistic jets to Lbol of their respective source. They con¬ 
sider both black-hole binaries (BHBs) and AGN which means that 
their study samples an extremely broad spectrum of both black 
hole mass [0.8 < log{M^}^/M q) < 10] and bolometric luminosity 
[37 < logLbol < 47]. Even though the jKing et al.H2013] l study deals 
primarily with the low-velocity soft X-ray warm absorber they do 
include a few ‘ultra-fast outflows’ (which they define as those with 
Vw ^ 0.01 c, as we have done in this work) as a matter of compar¬ 
ison, although they do not fit them directly as part of their regres¬ 
sion analysis; two of their ‘ultra-fast’ sources, namely 3C 111 and 


APM 08279-1-5255, are also included in our Suzaku sample. In prin¬ 
ciple we could compare the properties of the Suzaku-detected out¬ 
flows with those of King et al. to determine if our high velocity Fe K 
winds scale in a similar manner to the warm-absorber. However, 
this is not possible in practice because the two studies compute the 
wind kinetic power according to different assumptions (King et al. 
give their quantities in units of per-covering-fraction), which means 
the intercept of any linear regressions will therefore by intrinsically 
offset. Even so, the slope of any linear fit — which is the important 
parameter as it measures how the kinetic power varies with bolo¬ 
metric luminosity — can still be compared because it is largely 
independent of the normalisation. 

King et al. find the kinetic power of warm absorber to be 
strongly correlated with Lbol. scaling with a tight global slope of 
logLw (1.58 ± 0.07) logLboh while they found that those in AGN 
have a much flatter local slope of 0.63 ± 0.30 when considered 
in isolation (although that this could be driven by there being a 
smaller number of AGN). In Figure we showed that kinetic 
power of the Suzaku-detected FeK winds are also strongly cor¬ 
related with the bolometric luminosity with a slope in the range 
1.5^gg (Rp is 0.54). This result, which is driven by the v„ oc 
dependence, is in good agreement with the slope reported by King 
et al. for their global sample and suggests that the fast Fe K winds 
may scale with LeoI tti a similar manner to the warm absorber. To 
summarise, by comparing our results for those obtained by |King| 
jet al. 12013] we see that the fast Fe K winds (v* > 0.01 c) is similar 
to that found for slow warm absorbers. 


4.2 Are the winds radiatively accelerated? 


In this work we shown that the observed wind velocity, and subse¬ 
quently the overall wind energetics, are proportional to the AGN 
bolometric luminosity. It is therefore tempting to conclude that 
the winds are accelerated by radiation pressure. In this section we 
investigate whether radiation pressure alone can account for the 
observed correlations. In contrast to the undoubtedly line-driven 
winds which are prevalent in the UV spectra of AGN, the high 
ionisation state of the Fe K winds suggests that line-driving is un¬ 
likely to be the dominant acceleration mechanism here. An al¬ 
ternative means of radiatively accelerating high ionisation gas is 
through Thomson/Compton scattering of the continuum X-ray pho¬ 
tons. While line-driving may play a role in large black holes ac¬ 
creting near the Eddington limit (Hagino et al. 2014, in prep) we 
concentrate on the latter mechanism here. 

In comparison to fast line-driven winds in the UV, the com¬ 
paratively small interaction cross-section in highly ionised material 
means that high-ionisation winds, such as those considered here, 
couple much less efficiently with the incident radiation and, as a 
result, require either much larger source luminosities or higher col¬ 
umn densities to achieve equivalent outflow velocities. A general 
characteristic of these ‘continuum-driven winds’ is that they are 
accelerated by photons scattering off free electrons in the absorb¬ 
ing gas which transfers a portion of the photon momentum to the 
material, hence causing a wind. This scenario has been considered 
extensively in the literature (e.g., King^^ounds|2CiO^ Eng|2005| 


|2010[ |Reynold^|2012[ [Costa, Sijacki & Haehnelt||2014^ . out¬ 

line the general theory again here in an effort to search for relevant 
relationships with Lboi which we can relate to the results of our 
regression analysis. 

The efficiency through which the photon momentum of the in¬ 
cident radiation field is transferred into the wind depends upon the 
plasmas optical depth to electron scattering, T, which determines 
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Figure 3. Scatter plots showing logLboi against the raw wind parameters: (a) log(r/cm), (b) log(Mw/gs ^), (c) log(Lw/ergs *) and (d) log(pw/gcms ^). 
Data points and plot components have the same meaning as in previous figures. 



Figure 4. More scatter plots, this time showing logLboi against the normalised wind parameters: (a) log(r/rs), (b) log(Mvv/^Edd)> (c) log(Lw/^Edd) and (d) 
log(pw/ pEdd)- Data points and plot components have the same meaning as in previous figures. 
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how many times each electron is likely to scatter in the flow. The to¬ 
tal electron scattering optical depth of a homogeneous wind viewed 
from infinity down to radius r is given by 


T = Ny{ CTp — Oj / n- dr = 


j: 




ATtbm^Vvjr' 


( 8 ) 


where (Tj is the Thomson cross-section for electron-scattering, 
n = Mv,/AKbm^Vvid is the electron number density in a mass- 
conserving spherical flow (c.f., §2.2^ and is a geometrical factor 
which takes into account the global covering of the gas. Normal¬ 
ising to the Eddington mass-accretion rate, Mndd = ^Edd/flc^ = 
OjTjc, then yields 


(9) 

fcflcvw R Mfidd’ 

as has been shown by |King & Pounds| ( |20Q3| ) (see also |King|2010| 
for a recent review). If we further assume that the measured flow 
velocity roughly corresponds to the escape velocity at the radius 
where the wind was launched, i.e., R = Rg^c = 2GMbh/v^, we have 


1 Vw 

IbT] C M-Edf 


( 10 ) 


which, after substituting for Msdd = ^^Edd/flc^ and re-arranging, 
becomes 


Pw 


= 2bx^ = 2bl^ 

C Ac 


( 11 ) 


where the geometric factor b = and we have 

substituted from A = Tbol/AEdd in ths final step to express the 
relation in terms of observable quantities. Equation (Hj there¬ 
fore shows that a wind accelerated by electron scattering should 
have an outward momentum-rate which is proportional to that of 
the incident radiation field. This result (and the associated deriva¬ 
tion) is directly equivalent to the one presented by |King & Pound^ 
( |2003^ . Note that while T has thus far been described as the op¬ 
tical depth to electron scattering it can in principal also account 
for additional sources of opacity in the flow, e.g., though bound- 
free or bound-bound absorption, which can also be further boosted 
by multiple scattering events. A reasonable way of testing for the 
continuum-driven wind scenario is to therefore see whether and 
Pboi = Abol/c follow a linear trend. Eigure|^ shows that the two 
are clearly correlated, with a slope of j3 = 1.2(^)0^. This is consis¬ 
tent with the order of unity expected from equation (|TTJ, and is also 
agreement with the slope of jSxmm ~ 1.6 ± 1.1 found by |Tomb^ 
1 2013^. Eigure|^ is a constraint diagram for log(/>„/pbol) 


et al. 


showing the deviations from the expected ratio of unity; in total, 
most of the sources (13/20; 65% of the sample) are formally con¬ 
sistent with the ratio of unity expected for a continuum-driving sce¬ 
nario. 

Another useful relationship can be obtained by multiplying 
equation l[TTJ again by Vw: 


, _ r _ j, ^ , 

L-w — OX ^Edd — ^ A ^bol? 

C AC 


( 12 ) 


i.e., the wind kinetic power is predicted to be proportional to the 
bolometric luminosity of the AGN. In Eigurej^ we showed that Lw 
and Lboi are indeed correlated (Table [^, although the slope (jS ~ 
1.5^Qg) appears to be slightly steeper than (but still marginally 
consistent with) that predicted from equation The steeper 
slope here is probably due to the Vw ~ relationship observed 
before. Indeed, taking this into account yields oc x Lboi = 
which is exactly what we observe in these data. 


Even so, it is worth noting that other acceleration mecha¬ 
nisms are not conclusively ruled out for these winds. If the winds 
are launched with zero initial radial velocity, then the terminal 
velocity of the wind in terms of the Eddington limit is: V[„f = 
^2GMBHMaunch(/Tf-bol/AEdd ^ 1), where vm is the local ter¬ 
minal velocity and /x is the force multiplier (= 1 for Thomson 
scattering). As a result, any wind accelerated purely by Thomson- 
scattering will not be able to greatly exceed the local terminal ve¬ 
locity without either: (a) additional sources of opacity in the flow 
(such that /x > 1), or (b) additional mechanisms transferring mo¬ 
menta to the flow. In this case, additional opacity in the flow may 
be attributed to the outflowing gas if it is stratified and contains 
clumps of lower-ionisation gas. This extra opacity could perhaps 
be associated with the partially-covering gas required by many of 
our models in Paper 1. Alternatively, magnetic processes may play 
a role in the initial acceleration which provides additional momen¬ 
tum to the flow, thereby allowing it to escape the system. Indeed, 
have been proposed along these lines in the literature (e.g., |Ohsuga | 
|et al.|200^|0hsuga|2012[|0hsuga & Mineshige|2014| It is not pos¬ 
sible to distinguish between these possibilities on the basis of ex¬ 
isting data. 


4.2.1 Estimating the wind opacity 


Considering again the continuum-driven scenario, an interesting 
thing of note from Eigurej^ is that while all of the fast Ee K winds 
are consistent with the unity ratio within the errors, all of the slow 
winds fall below the limit (even when only the upper-limits are 
considered). This possibly reflects subtle changes in the x/X ra¬ 
tio which allows for the possibility of continuum-driven winds in 
sub-Eddginton AGN, whilst also allowing for subtle differences 
in the coupling efficiency of the gas between the two groups. As 
mentioned in SjT] the AGN considered here are relatively tightly 
clustered around a mean A ~ 0.1 which means that we are unable 
to reliably test how the Eddington ratio affects the Pw/Pbol ratio. 
Nevertheless, we can still try to gauge how x changes across the 
sample. Erom equation we can see that when normalised to 
Eddington Pw/PEdd = 2bx, with b ~ Sl/Xn ~ 0.4 (see S2.2i and 


therefore that 2/> ~ 1 to within an order of magnitude. This means 
that p„/PEdd — "T, and that the Pw/PEdd ratio can be in principal be 
used as a proxy for the gas opacity when the wind was launched. 
Figure|^shows how log(pw/PEdd) varies with Vw The higher ve¬ 
locity systems are largely consistent with p* = pedd ('r = 1) within 
the errors, with a mean of (T)fast ~ 0.3, which is consistent with 
what is argued by |King & Pounds] ( |2003^ . In contrast, the slow sys¬ 
tems (blue) all appear to have Pw/z^Edd ratios around an order of 
magnitude lower. Taken at face value, this possibly indicates that 
they have a lower opacity to scattering. There also appears to be a 
transitionary phase between the two groups, with Vw ~ a few thou¬ 
sand km *, which may be consistent to the idea of stratified winds 
launched over a wide range of radii with varying velocities (see 
|Tombesi et al.|20T3) ). It is worth noting, however, that the opacity 
estimated from Pw/PEdd could be lower than that of the bulk flow if 
(a) we view the wind from an acute angle where the observed veloc¬ 
ity is not representative of the bulk flow velocity, and (b) the wind 
has a component of completely ionised gas which contributes to the 
scattering opacity but has no spectroscopic signatures. One possi¬ 
bility is that all of the winds have T ~ 1 when integrated over all 
lines of sight and across all ionisation phases. Alternatively, the fact 
that the slow winds have T ^ 1, coupled with the fact that they are 
preferentially observed in lower luminosity AGN, may indicate that 
they are accelerated by another mechanism such as magnetic pres- 
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Figure 5. Comparison of the outflow (p„) and the radiation field (pbol) 
momentum rates. Top: a scatter plot showing the clear linear relationship 
between the two quantities (see Table|^for parameters). Plot components 
are the same as previous. Bottom'. Constraint diagram showing the p„/pboi 
ratio for each source. Almost all of the AGN in the sample are consistent 
with a ratio of unity (denoted by the dotted black line). 



log(v.j^/kln s ^) 


Figure 6. Plot showing log{pw/PEdd) versus log(vw/kms“*) which can 
be used as a diagnostic for T (see text for details). The majority of fast 
AGN (v„ > 3000 km s“'; red points) are consistent with pw = PEdd , thereby 
implying that T = 1, while all of the slow winds (blue points) are well below 
a ratio of unity. 


sure, for example, which may have a larger contribution at lower 
luminosities (e.g., see |Ohsuga et al.|2009[|0hsuga|2012^ 

Regardless, the fast winds — which are the ones most likely 
to affect their host galaxy due to their large velocity and therefore 
the most interesting in terms of AGN-host-galaxy feedback — are 
consistent with a wind launched by continuum scattering at the Ed¬ 
dington limit (e.g., [King & Pounds|2003^. 


4.3 Energetic significance and feedback implications 

The independent Suzaku and XMM-Newton wind samples have 
conclusively shown highly ionised Fe K winds are both frequently 
observed in AGN and that they contain a large amount of mechan¬ 
ical energy. However, a fundamental question still remains: is the 



Figure 7. Diagnostic plot comparing the ratio of L*. /L(,oi to the Lboi of each 
AGN in the Suzaku sample. The horizontal lines denote L„/Lboi fractions of 
0.5% and 5% which are the theoretical thresholds for feedback. Almost all 
of the fast absorption systems (red points) have L„ in excess of 0.5% of Lboi 
which suggests that they may be significant in terms of AGN-host-galaxy 
feedback scenarios. 


energy imparted by the mass flow likely to play a role in terms 
of AGN-host-galaxy feedback scenarios? AGN feedback models 
have postulated that ~ 5% of an AGNs’ bolometric radiative out¬ 
put needs to be converted to mechanical power in order to have a 
notable effect on the growth of a central SMBH and the host galaxy 
(e.g., [Di Matteo, Springel & Hem quist|2005 l. However, recent nu¬ 
merical simulations by Hopkins & Elvis| 2010| have shown that 
should the effects of ‘secondary’ feedback such as cloud ablation be 
taken into account, the required energy can be significantly lower, 
of the order ~ 0.5%. Moreover, hydrodynamic simulations have 
also shown that even relatively modest FeK outflows (e.g., those 
with only Vw ~ 0.01 c, Mw ~ 0.01 (Hq yr^* andLw ~ 10^^ erg s^*; 
see jWagner, Umemura & Bicknell||20f^ can impart significant 
feedback upon the host galaxy on ~kpc scales should these ef¬ 
fects be taken into account. Investigating what fraction of an AGNs’ 
bolometric output is conveyed through mechanical power therefore 
provides a means of qualitatively assessing the likely energetic sig¬ 
nificance of the mass flow. A plot of Lboi ^ (^w/Lbol) to this effect 
is shown in Figure]^ 

A large fraction of the outflows have a kinetic luminosity ei¬ 
ther in excess of or consistent with the 0.5% of Lboi, with 10/20 also 
exceeding 5% as well. The mean kinetic power considering only 
the fast winds is (L„/Lboi) ~ 7%. This is well above the supposed 
~ 0.5% threshold for significant feedback. In the most conserva¬ 
tive case, i.e., by considering only L™" given the available errors, 
we estimate (L^/Lboi) ~ 0.2%. This is comparable to the ~ 0.5% 
threshold and consistent with what was found by |Tombesi et al.| 
( |2012a| l on the basis of the XMM-Newton outflow sample. There¬ 
fore, even in the most conservative case, the energetic output of 
the fast Suzaku-delected Fe K winds is comparable to the threshold 
thought necessary for feedback. Importantly, these estimates are 
derived solely from the wind kinematics and are therefore indepen¬ 
dent of the initial driving mechanism. 


4.4 Relation to Mbh cr* in quiescent galaxies 

The FeK winds detected here may therefore impart enough energy 
into the host galaxy to play a roll in galaxy-scale feedback. With 
this in mind, it is useful to highlight the recent work of |McQuillin | 
|& McLaughlin|j2013| >. To set the context of their work, it is first 
useful to explain the difference between momentum- and energy- 
driven outflows. As a wind propagates into the host galaxy it will 
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Figure 8. Histogram comparing the Suzaku-observed distribution of out¬ 
flow velocities (blue cross-hatched area) to those inferred from local qui¬ 
escent galaxies by |McQuillin & M cLaughlm|^2013| who assume that their 
observed Mbh “ (7» relation (c.f., |Gu!teEi^^L|2QQ^ is induced by energy- 
driven feedback, with Mbh^’w C7^. The blue and black arrows correspond 
to the median outflow velocites of Vw ~ 0.056 c and 0.035 c in the Suzaku- 
observed AGN and quiescent sub-samples, respectively. The two samples 
are formally indistinguishable by the Kolmogorov-Smirnov (KS) test. See 
text for further discussion. 


sweep up a shell of ambient gas, the dynamics of which are de¬ 
termined by how efficiently the shocked wind material behind the 
shell can cool. Initially, the shocked gas can cool efficiently (by 
Compton scattering) so the shell stays thin and is driven outwards 
directly by the force of the wind (cf. our equation pTj also |King| 
|& Pounds| ( |200^ ; [KingH20l^ ; |McQuillin & McLau^lin| ( |2012^ ) 
This is momentum-driven feedback. As the outflow reaches larger 
distances, however, the cooling efficiency drops and the shocked 
gas becomes geometrically thick and hot. The outwards force on 
the shell then comes from the wind thrust in equation (U} medi¬ 
ated by the thermal ppressure and work done on and by the shocked 
gas. This is energy-driven feedback (e.g., |King|2003| 12005] |Zubo-| 
|vas & King|2012[[McQuillin & McLaughlin|2013^ . The transition 

between the two phases can occur at relatively small radii (see|Mc-| 
IQuillin & McLaughlin|2013] , and energy-driven outflows are there¬ 
fore expected to be dominant on large galactic scales. 

IMcQuillin & McLaughlin|2013| consider purely energy-driven 
outflows and investigate their link to the observed Mbh ^ cr* rela¬ 
tionship. They showed that, for an energy-driven shell to reach the 
escape speed from an isothermal sphere of dark matter and gas with 
velocity dispersion a requires a critical BH mass given by 


/VwN 


1 c; 

(lOSMey 


= 6.68x 10^^ 


200 km s^* 


5 


(13) 


such that the velocity of the wind from the black hole (vw) enters 
explicitly. McQuillin & McLaughlin used this equation to infer the 
wind speed that would have been needed during the active proto- 
galactic phase of 51 quiescent {z ~ 0) galaxy spheroids (taken from 
the sample of |Gultekin et al.|p009| >) in order to explain the ob¬ 
served values of Mbh ™d a in these galaxies as a result of (energy- 
driven) feedback. Their main result is a distribution of BH wind 
speeds, derived in essence from the scatter of the Mbh ^ ela¬ 
tion, that bears a remarkable resemblance to the empirical distribu¬ 
tions of Vw in the AGN samples of both |Tombesi et al.|2010a| and 
Paper I. 

Figure [8| compares the velocity distribution inferred by |Mc-| 
IQuillin & McLaughlinI ( |2013[ ) to the one that we measured for 
SMZflfcn-detected Fe K outflows (note that we consider all of the out¬ 
flows in this panel). As noted by McQuillin & McLaughlin, there 


is a close resemblance between the two, with both samples tracing 
a similar overall distribution. Their overall range are similar, as is 
their median: Vw ~ 0.035 c in the quiescent galaxies, and ~ 0.056 c 
for the Suzaku-detected outflows (see Paper I). Formally, the two 
distributions cannot be distinguished: a Kolmogorov-Smirnov (KS) 
test only rules their being from the same parent population at 
Pks ~ t5% confidence. The fact that there are similarities between 
the inferred and observed Vw distributions at all is in itself ex¬ 
tremely interesting, not only because the independent samples of 
quiescent galaxies and AGN are drawn from two different sub¬ 
populations of galaxy, but because it implies a close relationship 
between the emergent wind and the evolution of the host galaxy. 
Overall, the work of |McQuillin & McLaughlin|]2013^ is encourag¬ 
ing for the general notion that high-velocity outflows contribute to 
the feedback which is invoked to explain the Mbh ^ relationship 
in local quiescent galaxies, and, more importantly, for the idea that 
AGN play a key role in galaxy formation and evolution. 


5 SUMMARY & CONCLUSIONS 

In Paper I we formed a detailed spectroscopic study of Fe K ab- 
soption in a heterogeneous sample of 52 Suzaku-ohserved AGN. 
There, we found that ~ 40% of the sources in the sample harboured 
statistically significant Fe XXVI Hea and/or Fe XXVI Lya absorp¬ 
tion lines with velocities ranging from a few thousand kms^* to 
~ 0.3 c. We also measured their column density ionisation parame¬ 
ter of the absorber using the XSTAR photoionisation code. Here, we 
build upon these results to compute the location r, mass outflow rate 
Mw. kinetic power Lw and momentum flux pw of the implied wind, 
and use these results to inform an exploratory analysis of how the 
wind parameters scale with the bolometric luminosity of their host 
AGN. We focus our attention on the fastest winds, defined here as 
Vw > 0.01 c, because these ones likely represent genuine disk-winds 
which are likely distinct from the slower (and more distant) warm 
absorber. The main results of our study are summarised below. 

(i) The mean radial distance to the fast winds is 
(r) ~ 10^^^*^cm (typically ~ 10^^^rQ from the black 
hole. This corresponds to ~ 0.0003 — 0.03 pc, such that 
the FeK absorbers are located much closer to the black 
hole than the traditional parsec-scale warm absorber. The 
mean wind mass outflow rate and mean kinetic power are 
constrained to (Mw) ~ 10^^^^® gs^* (~ 0.01 — IMwyr^^) 
and (Lw) ~ 1043 -44 

erg s *, respectively. The average upper 
limit is ~ MEdd> but is still an appreciable ~ 1% of Medd 
if only the lower limits are considered, whereas the kinetic 
power is constrained to ~ (0.01 — 0. l)LE(jd- These properties 
are consistent with wind at a substantial fraction of the 
Eddington limit which carries a large amount of mass into 
the host galaxy. 

(ii) The Fe K absorber column density and ionisation parameter 
are largely independent of Lbol but there is a significant 
correlation between the wind velocity Vw and Lbol> with 
a slope P = 0.4+Q2. This tentatively indicates that more 
luminous AGN launch winds with a larger observed velocity 
along the line-of-sight. 

(iii) The winds in more luminous AGN contain more material, 
and they are subsequently more energetic. The median slope 
of the L(,ol ^ t-v/ and L(,ol ~ Pw regressions are formally 
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consistent with unity. 


(iv) There are no significant correlations in the L(,ol ^ (''/^s) and 
^bol “ (j'^w/^Edd) planes. Therefore, the winds are located 
at similar distances in Schwarzschild units and contain a 
similar fraction of Medd- 


(v) The momentum flux of the radiation field (pbol) and that 
of the wind (p„) are strongly correlated, with a slope 
jS = 1.2^g'^. This is quantitatively consistent with what is 
expected from a continuum-driven wind. 


(vi) A significant fraction of the sample (17/20, 85%) exceed the 
minimum Lw/Tbol ~ 0.5% threshold thought necessary for 
feedback (Hopkins & Elvis|2010 1 , while 9/20 (45%) also ex¬ 
ceed the less conservative ~ 5% threshold as well jPi Mat-| 
|teo, Springel & Hernquist||2005[ l. In the most conservative 
case, the mean L„/Lbol ratio is ~ 0.2%. This suggests that 
the winds may be sufficiently energetic in terms of feedback. 


These results enforce those recently obtained with XMM-Newton 
and provide additional evidence in favour of Fe K absorption being 
a prevalent feature of the AGN X-ray spectrum, and that the ensu¬ 
ing wind may be energetically significant in terms of feedback. 
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